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■ We analyze the Blume-Emery-Griffiths model with disordered magnetic interaction displaying the inverse freezing phenomenon. The 
. ^ behaviour of this spin-1 model in crystal field is studied throughout the phase diagram and the transition and spinodal lines for the 

O model are computed using the Full Replica Symmetry Breaking Ansatz that always yelds a thermodynamically stable phase. We compare 

_ the results both with the quenched disordered model with Ising spins on lattice gas - where no reentrance takes place - and with the 

* I , model with generalized spin variables recently introduced by Schupper and Shnerb [Phys. Rev. Lett. 93, 037202 (2004)]. The simplest 

$_( , version of all these models, known as Ghatak-Sherrington model, turns out to hold all the general features characterizing an inverse 

'^T^ ■ transition to an amorphous phase, including the right thermodynamic behavior. 

B- 
^: 

B 1 Introduction 
I 

In the last years, in a relatively broad number of materials an apparently weird kind of phase transition 
O has been detected and investigated: the "inverse melting". Such a phenomenon was, actually, already 
hypothesized by Tammann a century ago [J but it has now been found experimentally in polymeric and 
colloidal compounds, high-Tc superconductors, proteins, ultra-thin films, liquid crystals and metallic alloys 
|Si3, 4, 5, 6, 7, 8^. This kind of transition, includes, e.g., the solidification of a liquid or the transformation 
• of an amorphous solid into a crystal upon heating. The reason for this counter intuitive process is that a 
t phase usually at higher entropic content happens to exist in very peculiar patterns such that its entropy 
[ is decreased below the entropy of the phase usually considered the most ordered one. An example taking 

■ place in the widely studied Crystalline Poly(4-methylpentene-l) (P4MP1) |5] is the one of a crystal state of 
[ higher entropy that can be transformed into a fluid phase of lower entropy on cooling. Inverse transitions, 

in their most generic meaning (i.e. both thermodynamic or occurring by means of kinetic arrest), have 
been observed between fluid and crystal phases 0, between glass and crystal j2j and between fluid and 
glass ("inverse freezing") |^. For a comprehensive summary of materials undergoing some kind of phase 
transition one can turn to Ref. [Hj. The aim of this work is to study a simple mean- field model for the 
inverse transition in a spin-glass, in order to heuristically represent the inverse fluid-amorphous transition. 



. ^ 2 The Random Blume-Emery-GrifRths-Capel Model 

X; 

. We have been analyzing the Blume-Emery-Griflfiths-Capel (BEGC) model with quenched disorder using 
the Full Replica Symmetry Breaking (FRSB) scheme of computation that yields the exact stable thermo- 
dynamics |1H I12 | [T3]. It includes the Blume-Capel 14, and the Blume-Emery-Griffiths |15) models, when 
the couplings are ferromagnetic, and the Ghatak-Sherrington (GS) model, when the couplings are 
random variables and no biquadratic interaction occurs. 
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The model Hamiltonian we consider is 

^ K 

ij i=l i<j 

where S" = 1, 0, —1, D is the crystal field, Jij are quenched random variables (Gaussian) of mean zero and 
variance 1/N. The parameter K represents the strength of the biquadratic interaction. A generalized spin 
variable has been recently proposed jlUl lH]. for which the degeneracy of the magnetically interacting sites 
{S = ±1) can be larger (or smaller) than the one of the holes (S = 0). To be as inclusive as possible we 
have been looking at the phase diagram of the random BEGC model also in terms of these generalized 
variables considering the thermodynamically stable spin-glass obtained by means of the FRSB Ansatz. 

Let us call call k the degeneracy of the filled in sites of one type (5 = 1 or 5 = —1) and / the degeneracy 
of the empty sites {S = 0). The relevant parameter is r = k/l. When r = 1 the spin-1 model is obtained. 
If, furthermore, K = 0, the model is the GS one. When, otherwise, r = 1/2 and Z? — /i = —D the lattice 
gas formulation of Ref. is recovered, for which no reentrance was observed jlllll2j . 

We discuss the physically stable solution for both the GS model |161ll8j {K = 0) and the model with 
attractive biquadratic interaction jl2j (K/J = 1), whose phase diagrams are plotted in Fig. ^ and we 
compare it with the Replica Symmetric (RS) results (see Fig. [JJ. We study the behaviour of the phase 
diagram for (a) the lattice gas case (r = 1/2), for which no inverse transition occurs anywhere in the 
parameter space; (b) the spin-1 case (r = 1), where the first order transition line displays a reentrance 
soon below the tricritical temperature; (c) the generalized cases as r > 1, in particular we plot the results 
of the model with variables taking values 5* = {1,1,0, —1,-1} for which the reentrance takes place above 
the tricritical point, along the second order phase transition line. 



3 The Full Replica Symmetry Breaking scheme of computation for the Random BEGC model 

Using the replica trick [201 the free energy of n replicas of the system with Hamiltonian Eq. turns out 
to be 



n/3/ = § + I) X; + ^ E - log E -M-Pn') (2) 

^ ^ a=l a^b {S} 

^ ^ a=l a=l a^b 



(3) 

where the order parameters introduced in the computation satisfy the self-consistency equations 

Pa = {Sl) qab = {SaSb) (4) 

Choosing the Parisi Ansatz for the matrix qab and performing infinite steps of RSB ^H] the free energy 
of the spin-glass phase appears to be 

Pf = ^(K+^)p'-^ [' dx q\x) - Pm 0) (5) 



2V 27' 4^0 
and the self-consistency equations become 

/■°° , r./ X.// n2 /"°° , r... N 2re®cosh/3y 

qix) = J dy Pix,y)(f>'{x,yf p = J dy P{l,y)-^-^-^-^—^ (6) 
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Figure 1. The D-T phase diagram in absence of biquadratic interaction. Three models with different behaviors are plotted: r = 1/2, 1, 
2. For each r-model three curves are represented, each departing from the same tricritical point: the full curve on the left is the spinodal 
of the PM phase, the dashed one in the middle is the first order transition line and the right one is the SG spinodal line. The group of 

three curves on the left are for the Ising spins on lattice gas (r = 1/2, Tc = 1/3, Dc = 0.73105). The group of curves in the middle 
represent the lines of the GS model (r = 1, = 0.96210). The curves on the right correspond to the r = 2 model (Dc = 1.19315). In 
the inset the same diagram is plotted when K = 1 (attractive biquadratic interaction). No qualitative difference occurs. 

Figure 2. D-T diagrams a.t K = 0: Ising spin glass on lattice gas (left), GS model (center) and random BEGC model with r = 2 
(right). Both the RS and the FRSB solution (the latter with error bars) are plotted. For r = 1/2 no reentrance occurs. A reentrance 
occurs, instead, below the tricritical point in the GS model along the first order transition line (and a second reentrance seems is there 
for T < 0.03). In the latter model the reentrance occurs above the tricritical point. 

where Q = /?^/2[/3 — q{l)] + /3Kp — f3D. The functions (p{x, y) and y) are solutions of the non-hnear 
partial differential equations 

kx.y) = -| [</.'(x,2/)]'} (7) 

P{x,y) = -^[P"{x,y) -2(ix[P{x,y)(^'{x,y)\'] (8) 

with boundary conditions 

<A(l,y) = ^log(2 + 4re®cosh/3y) P(0,y)=%) (9) 

The dot and the apex respectively represent the derivative with respect to x and y. For a detailed study 
of the above equations the reader can look at Ref. I12j . 

Having a model with variables displaying a relative degeneracy r [D = D,.), in order to describe the 
partition function of another model whose variables have degeneracy r' it is enough to vary the crystal 
field as Dr' +Tlogr' = Dr + Tlogr. This does not hold, however, for the state functions obtained deriving 
the thermodynamic potential with respect to the temperature (e.g. entropy and internal energy) that will, 
instead, receive contributions from additional terms. For instance, the entropy function (K = 0) is 

s{P, D) = -^[p- q{l)f -^p[p- q{l)] - PpD + plog{2r) + dx q\x) (10) 

/CO 
dy P(l, y) log {2 + 4 exp (0 + log r) cosh fiy} 
-oo 



Identifying Z)i/2 = ~P i^i Eq (fTU]) one recovers the case of spins on a lattice gas of chemical potential 
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4 Thermodynamics 

By solving Eqs. (|5|)- (|lfl|l we build the D-T phase diagrams that we plot for K = {) and i^T = 1 in Fig. ^ 
and inset, at different values of r. The reentrance in the D-T plane is present already in the spin-1 GS 
model. As a consequence this implies that there is no need for the intuition of Ref. |inj in order to have 
a model for inverse freezing from low temperature liquid to high temperature amorphous solid. This is 
at difference with the liquid-crystal inverse transition ( "inverse melting" ) for the description of which the 
original Blume-Capel model is not adequate and r > 1 is needed |in[ l9l. 

The transition lines are not very much dependent on the Ansatz used to compute the quenched average 
of the free energy. For not extremely low T, the first order transition lines yielded by the RS and the FRSB 
Ansatz even coincide down to the precision of our numerical evaluation of the FRSB antiparabolic Parisi 
equation Q whereas they slightly differ at very low temperature (see inset of Fig. For what concerns 
the spinodal lines, the RS ones are shifted inside the pure paramagnetic (PM) phase. 

Introducing a biquadratic interaction term and varying it from attractive to repulsive the situation does 
not change much (see e.g. inset of Fig. For any value of K no reentrance of the phase transition line 
occurs in the D-T phase diagram of the lattice gas model whereas it is always there for the spin-1 model. 
The only consequence of reducing K is that the area of the phase coexistence region is reduced (the 
tricritical temperature tends to zero as, K ^ — oo). 

The slope of a first order line is given by the Clausius-Clapeyron equation. For the BEGC model it can 
be written in terms of the crystal field D (playing the role of a chemical potential), instead of the pressure 
that is not defined in our model:^ 



This formula is valid for any r. We stress that in passing from r to r' also the entropy (Eg. Ilflj) changes of 
a term p log r/r', in agreement with the crystal field shift given above. 

In Fig. 0the behaviour of the entropy as a function of the temperature (Eq. (|10|) ) is shown across an 
inverse transition (as a function of the crystal field D in the inset) for the spin-1 model. The entropy of the 
PM phase below the first order transition line is smaller than the entropy of the SG: heating the system 
the paramagnet becomes an amorphous magnet (i.e. "freezes" ) acquiring latent heat from the heat bath. 
Going down along the transition line, as As changes sign in Eq. Hll|) the slope becomes positive. The 
As = point is called a Kauzmann locus 

A reentrance in the transition line can be due both to the existence of a liquid phase with an entropy 
lower than the one of the solid or to a the liquid phase more dense than the solid one (like in the water-ice 
transition). When an entropic inversion accounts for the phase transition, the equilibrium transition line 
changes slope in a point where the entropy of the fiuid phase, suq, becomes equal to the one of the solid, 
Ssoi, according to the Clausius Clapeyron equation for first order phase transitions. From this point on 
a whole iso-entropic line. As = suq — Sgoi = 0, can be continued both inside the solid and the liquid 
phases. This is a particularly interesting observation since, in the context of glass formers, a transition to 
an "ideal" glass at the temperature at which As = is hypothesized, in order to avoid the paradox that 
an under-cooled liquid might possess less entropy than the associated crystal at the same values of the 
thermodynamic parameters (21] . From an experimental point of view this Kauzmann temperature would 
be the temperature of the glass transition (that is not a true phase transition because strictly kinetic in 
origin) in an idealized adiabatic cooling procedure. Since the astronomically long relaxation time needed 
to actually perform such an experiment makes such a procedure unfeasible, the evidence in favour of the 
existence of a thermodynamic glass transition mainly comes from analytical and numerical investigations 
(see e.g. Refs. j232S|)- The fact that a As = line turns out naturally in the description of the behavior of 
materials with inverse transition avoids, at least for these substances, the Kauzmann paradox and breaks 
the connection between As = extrapolation and the existence of an ideal amorphous phase |31I21- 




dT ppm - Psg Ap 



^With respect to the classic Clausius-Clapeyron equation D takes the place of the pressure and p plays the role of the specific volume. 
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Figure 3. Free energy -Fpm and F^g^^ vs. D at T = 0.23 in the GS model (K = 0). The left and right side vertical lines are at the 
spinodal points of, respectively, the PM and the SG phase. The first order transition occurs at Dj^at = 0.9344. In the inset 

F^^ is displayed. At this temperature the two functions merge very near to the tricritical point (right side vertical line). 

Entropy vs. T at the crystal field value of _D = 0.96 for the GS model. Tjst = 0.302. The two vertical dashed lines are the SG 
(TsG = 0.314) and PM (Tpm = 0.3465) spinodal lines. In the inset s{D) is plotted at T = 0.23. 
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5 Conclusions 



We have shown that the Ghatak-Sherrington model undergoes the inverse freezing phenomenon acquiring 
latent heat from the heat bath as the paramagnet becomes a spin-glass. Many other models can be built 
starting from this one, introducing an attractive or repulsive biquadratic interaction (the last term in 
Hamiltonian (^) and/or tuning the relative degeneracy of the value S = and S*^ = 1 of the spin 
variable (as in ^Ul) but the GS one already contains all the features needed to qualitatively represent the 
experimental results. 

We have considered here a case where the frozen phase is a spin-glass and the transition is only static: 
in order to have a structural glass and to deal with a dynamic (mode coupling) transition, it would be, 
however, enough to apply the simpler one step RSB Ansatz in, e.g., a model with a p-body magnetic 
interaction (with p > 2 instead of a couple interaction), such as the one introduced in Ref. i24| or an 
analytically solvable model with spherical variables instead of discrete ones and again a multi-body inter- 
action (a valuable starting point would be the work of Ref. [25\). Else, a spin-1 model with an orthogonal 
matrix Jij, instead of the Gaussian one, can be used, as latterly done by Sellitto [26j . Besides the inverse 
freezing phenomenon this last model also displays a fluid-fluid transition (between two distinct paramag- 
netic phases), a feature recently found in the solution of 4-methil-piridyne (4MP) and a-cyclodextrin in 
water, investigated by Angelini and coworkers jS|. 

With the help of the class of models presented here the connection between entropy driven phase reen- 
trance and shear thickening can also be tackled ;27' and, furthermore, a generalization of the spin-1 variable 
to a composition of "fast" and "slow" variables ^ coupled to two different thermal baths allows for studying 
anomalous latent heat in out of equilibrium transitions j28| . 
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